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Abstract. We review the notion of reducibility and we introduce and discuss 
the notion of orbital reducibility for autonomous ordinary differential equations 
of first order. The relation between (orbital) reducibility and (orbital) symmetry 
is investigated and employed to construct (orbitally) reducible systems. By 
standard identifications, the notions extend to non-autonomous ODEs of first and 
<^ ■ higher order. Moreover we thus obtain a generalization of the lambda symmetries 

of Muriel and Romero. Several examples are given. 
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1. Introduction and preliminaries 

In the present paper we discuss reducibility and reduction for ordinary differential 
equations. Our focus of interest is on (explicit) reducibility imparted by some map 
to a lower-dimensional system (possibly defined on a submanifold of n-space). For 
equations of higher order, reduction is frequently understood as reduction of order, 
which will also be considered. The main purpose of the paper is to consider various 
notions of reducibility, establish a general framework, and elucidate the relations 
between symmetry, generalizations such as A-symmetry, and reducibility. We 
clarify and extend notions, generalize results, and obtain new applications. In a 
related paper [2] we present a more thorough discussion of higher order equations, 
including prolongation formulas. 

For autonomous first order differential equations (resp. the associated vec- 
tor fields) there are well-defined and well-established notions of symmetry (sending 
parameterized solutions to parameterized solutions) and orbital symmetry (send- 
ing solution orbits to solution orbits, and consequently invariant sets to invariant 
sets). A canonical notion of reducibility (which includes symmetry reduction) was 
introduced and discussed in [3], and we extend this by introducing the notion 
of orbital reducibility. The latter turns out to correspond to a generalization of 
A-symmetries (Muriel and Romero; in particular [S]). Due to a relation between 
(orbital) reducibility and (orbital) symmetry there is a canonical construction of 
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(orbitally) reducible systems from (orbitally) symmetric ones, and in some cases 
it can be shown that all (orbitally) reducible systems are obtained in this man- 
ner. For vector fields the notions of " reducibility via some map" discussed in this 
paper seem to be comprehensive. Moreover, all notions of reducibility for non- 
autonomous systems or for equations of higher order can be traced back to the 
case of first order autonomous systems, just as non-autonomous systems or systems 
of higher order can be written as first-order autonomous systems. Reducibility for 
a non-autonomous first-order system (via a map sending solutions to solutions of 
some lower dimensional system) amounts to orbital reducibility of an associated 
autonomous system. This general framework for reduction comprises the main 
focus of the paper. 

Reducibility of a higher-order equation (in the sense of reducing order) is 
equivalent to reducibility of the canonically associated first-order system (in the 
sense of reducing dimension). In the course of proving this, we note that every 
(non-autonomous) m-dimensional system of order one with nontrivial right-hand 
side may be rewritten as a single equation of order m + 1, reverting the usual 
procedure. 

The problem to explicitly determine a reducing (or orbitally reducing) 
map for a given autonomous equation of first order seems to be just as hard 
(and as algorithmically inaccessible) as the problem of finding a symmetry (or 
orbital symmetry). This is essentially due to the straightening theorem and the 
implicit function theorem, which are not, or not completely, constructive. But the 
inverse problem to determine all vector fields reducible by a given map is easier to 
access in some relevant cases. In particular we transfer the approach from Olver 
and Rosenau [13] to the ordinary differential equation setting, and determine all 
differential equations which admit reduction by invariants of a given compact and 
connected group. 

Moreover we construct reducible higher-order equations from equations 
admitting symmetries, in particular Lie point symmetries. In this way we extend 
the class of reducible equations obtained from lambda-symmetric systems, and 
provide a different perspective for the latter. 

Throughout this paper we restrict attention to analytic functions and vec- 
tor fields; many of the results can, with some care, be extended to the smooth 
case. In order to give a self-contained discussion, and to make the paper acces- 
sible to readers with different backgrounds, we include a review (and sometimes 
rephrasing) of some facts and methods. 



We first fix some notation. Let an (analytic) autonomous ordinary differential 
equation 



be given on the open and connected subset U of IK" (with IK standing for M or 
C). We denote by Xj the corresponding Lie derivative which acts on analytic 
functions via 
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and recall that cf) is called a first integral of ([T]) if Xf[(f)) = 0. We distinguish 
between solutions (including parameterization) and solution orbits (trajectories). 
If one is primarily interested in orbits, then it is appropriate to consider an 
equivalence class of differential equations, rather than the single equation ([T]) in 
case / 7^ 0. Let Z be the zero set of /. Then a differential equation defined on 
the open set U ^ U has the same solution orbits on ?7 \ Z if and only if it has the 
form 

X = fi{x) ■ f{x) on U \ Z, 

with fi : U\Z ^ K analytic and without zeros. (The non-obvious direction holds 
because a suitable reparameterization of a solution of ([T]) will produce a solution 
oi X = fi{x) ■ f (x) .) Thus one has the notion of local orbital equivalence for vector 
fields that are defined on some open and dense subset of U . 

Remark 2.1. Two equations are locally orbit- equivalent if and only if they admit 
the same first integrals near any non- stationary point. Stated in a different way, 
two equations are locally orbit- equivalent if and only if they admit the same local 
invariant sets near any non- stationary point, o 

Remark 2.2. For a non- autonomous equation 

X = q{t, x) on V CKx K" 
one may define the " autonomized" system 

^° N , {xo,x) G V. 

X = q{Xo,x) V 5 / 

A reverse to this procedure is obtained as follows: If 

( /i(^) \ 

V fn{x) j 

and one of the components, say fi, has no zero on the open subset U , then passing 
to an orbit equation 

dx2/dxi = /2(x)//i(x) 

dXn/dXi = fn{x)/fl{x) 

provides a non- autonomous system whose autonomization is locally orbitally equiv- 
alent to ([1]) on an open and dense subset of U . Thus non- autonomous equations 
in dimension n and (local) orbital equivalence classes of autonomous equations in 
dimension n + 1 stand in correspondence, o 



2.1. Symmetries and orbital symmetries: Review. A symmetry of the 
autonomous differential equation ([1]) is a (locally invertible) map sending param- 
eterized solutions to parameterized solutions. An orbital symmetry of ([T]) is a (lo- 
cally invertible) transformation mapping solution orbits to solution orbits, hence 
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sending ([T]) to an orbit-equivalent equation x = fi{x) ■ f (x) . We recall a charac- 
terization of infinitesimal (orbital) symmetries; see e.g. Olver [12], Chapter 2, in 
particular Exercise 2.19, or see |16] . 

Proposition 2.3. Let g be a vector field on some open subset of U . Then: 

(a) The local transformation group generated by g consists of local symmetries of 
i = fix) if and only if [g, /] = . 

(b) The local transformation group generated by g consists of local orbital symme- 
tries of X = f{x) if and only if [g, f] = a ■ f for some scalar function a . 

For non-autonomous equations 

x = q{t,x) on\/CKxK" 

the usual definition of a symmetry is that of a locally invertible map defined on 
some open subset of V (thus transforming both t and x) such that solutions are 
mapped to solutions, see e.g. Olver [12]. Equivalently, by the above Proposition 
and Remarks 12.11 and 12. 2[ such a map is an orbital symmetry for any autono- 
mized system. Thus, necessary and sufficient conditions defining an infinitesimal 
symmetry of a non-autonomous first-order equation are known. 

2.2. Reducibility: Basic notions and results. The intent underlying any 
notion of "reducibility by some map" for equation ([T]) is rather obvious, as noted 
in the Introduction. But details have to be specified. The following particular 
form was introduced in [3]. 

Definition 2.4. We call the equation ([1]) reducible on U if there exist a number 
m, < m < n, a positive integer r, an analytic map \l/ : U — t- IC, and a 
differential equation 

y = h{y) 

defined on an open neighborhood U of "^{U) such that \& maps parameterized 
solutions of a; = f{x) to parameterized solutions of y = h{y) , and the derivative 
D\E'(x) has rank < m on f/, with rank = m at some point. 

The solution-preserving property is equivalent to the identity 

L>^(x)/(x) = /i(*(x)) on U. (2) 

Moreover, due to our assumptions the derivative D\E'(x) has maximal rank m on 
an open-dense subset of U . 

This definition of reducibility is designed to include interesting cases, like 
reduction by group invariants, which do not a priori provide a map to a vector 
space of smaller dimension. But locally near any maximal rank point, one has 
reduction to K™, and the structure of reducible vector fields is quite simple. In 
this sense, the problem can be locally "trivialized". 

Lemma 2.5. Let y ^ U such that D\l/(?/) has rank m, and let ipi, . . . ,ipr 
denote the entries of . Then there is a neighborhood U of y such that up to 
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a coordinate transformation one may assume that Tpi = Xi, . . . ,t(jrn = Xm, o-nd 
\I/ := {tpi, . . . , tpm)^'^ is then a reducing map to K"^ . Moreover, up to this coordinate 
transformation one has 



/ fi{Xi,. . . ,Xm) \ 
f m (-^1 ; • • • ) •^m) 



where the asterisks symbolize functions of all variables Xi, . . . ,Xn- 



Proof. We may assume by the implicit function theorem and the rank condition 
that ipi = Xi, . . . , 1pm = Xm. Since the matrix {Dipi, . . . , DiprY^ has rank m, the 
functions ipm+i, ■ ■ ■ y'lpr depend on xi, . . . ,Xm only. Then by the reducing property 



fi,...,fm can depend on xi, 



only, and all assertions follow. 



Remark 2.6. (a) This result can be refined. According e.g. to fTdf . Prop. 3.4 

every set 

{x e U; rank = q} 

is invariant for x = f{x) and locally a submanifold of U . Restriction to such sets 
thus suggests, and in principle allows, a case-by-case approach to reduction, 
(b) The problem of reducibility can be "trivialized" in yet another way: By the 
straightening theorem, f may be transformed to a constant vector field near any 
non- stationary point, for which reducibility (to any dimension > 1) is obvious. 



This observation shows that the set U in Definition 2.4 may play an important 
role. Furthermore, explicit determination of a reducing map for a given equation 
(or explicit determination of all vector fields reducible by a given map) is a different 
matter, and reducibility is a nontrivial property near stationary points (see also 

W- o 

The above trivialization results are of little practical relevance, due to their 
reliance on non-constructive theorems. Their principal value lies in providing 
insight into the local structure of reducing maps and reducible vector fields. 

We next rephrase and generalize some results of |3], Section 2, about the 
correspondence between reducing maps and involution systems. It seems appropri- 
ate to start with a relatively abstract statement, to clarify the relevant properties 
of the underlying algebraic structures. To motivate the role of function algebras 
in the following Theorem, note that for a given reducing map \1/ = (V'l, • • • j'ipr)^^ 
according to Definition 12. 4^ the algebra of all functions p{ipi, . . . , ipr) , P analytic in 
r variables, will be mapped to itself by Xf . Thus function algebras are naturally 
associated to reducing maps. 

Theorem 2.7. Given the analytic differential equation x = f{x) on U , let 
U U be open and connected, and A{U) the algebra of analytic functions from 
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U toK. 

(a) Let A4. he a Lie algebra of vector fields on U and denote by A{U)Ai the A{U)- 
module generated by M. . Let I{M.) C A{U) denote the algebra of invariants of 
A4. , thus Xg{(j)) = for all g ^ M. and all G I{.M.) ■ Then for any vector field 
f on U one has 

[/, M] C A{U)M =^ Xf{I{M)) C I{M). 
In particular, if f normalizes Ai then Xf stabilizes 

If, moreover, there are finitely many ai,...,ar € I{M) such that every element 
of I{M) can be expressed as an analytic function of the aj then one obtains a 
reducing map (ai, . . . , a^)*'^ for x = f{x) . 

(b) Conversely, let B be a subalgebra of A{U) , and f a vector field such that 
Xf{B) C B. Then f normalizes the Lie algebra 

C{B) = {g-Xg{B) = Q}. 
Moreover C{B) is a module over A{U) . 

If, moreover, this module is finitely generated, say by gi, . . . ,gs, then the gi are in 
involution on U , thus there are ^ij^ G A{U) such that for all i and j relations 

[9i,9j] = ^fJ^ijkgk 

k 

hold. 

Proof. To verify the nontrivial assertion of (a), let g E Ai and ip G I{Ai). 
Then by hypothesis, 

= X[,j^{^p) = X.Xjiij) - XfX.i^) = XgXji^) 

and therefore Xflip) G I{M.). Concerning (b), note that for every p E B and 
every g G C{B) one has 

X[,j]{p)=X,Xf{p)-XfX,{p) = 

in view of p E B, Xf{p) E B . ■ 

Remark 2.8. (a) There are obvious modifications of the Theorem for germs 
of local analytic functions and vector fields, resp. for polynomial and rational 
functions and vector fields. 

(b ) An important class of examples is formed by the systems symmetric with respect 
to a Lie algebra M. , thus [/, M] = 0. o 

Locally, the finite generation property holds in many cases, but such results 
are only partly constructive: 

Corollary 2.9. Let the analytic differential equation x = f{x) be given on U , 
and let y E U . 
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(a) // there is a rank s analytic involution system gi, . . . ,gs in a neighborhood of y 
(thus the gi are in involution and gi{y), . . . ,gs{y) span an s -dimensional subspace 
of K."- ), and there exist analytic functions Xij such that 



then there is a local reducing map \E', whose entries are common first integrals of 
the gi , to some equation on an open subset of K""** , with rank n — s . 
(b) // there is an analytic reducing map ^ as defined in ([2]), and D'^{y) has 
maximal rank m, then there is an analytic involution system of rank n — m, 
defined in some neighborhood of y , such that the entries of \1/ are common first 
integrals of this involution system, and the identities from part (a) hold. 

Proof. Part (a) is a direct consequence of Frobenius' theorem (see e.g. Olver 
[T2] . Section 1.3). For part (b) consider, in a suitable neighborhood of y, the 
homogeneous system of hnear equations 



(over the quotient field of the ring of analytic functions) and determine a basis 
(?!,... ,gn-m of the solution space. Obviously one may choose a basis consisting 
of analytic vector fields. Since the [giigj] are also solutions of this linear system, 
they are linear combinations of (?!,... , gn-m- The identities involving [g^, f] follow 
from Theorem 12. 7[ in view of the fact that any common first integral of the gi is 
locally a function of the ipj ; see Lemma 12.51 ■ 

Remark 2.10. The proof of part (b) shows that - in contrast to Frobenius - 
a corresponding involution system can be determined explicitly from the reducing 
map. o 

To finish this subsection, we discuss the relation between reducible systems 
and symmetric ones. 

Proposition 2.11. Let the analytic differential equation x = f{x) be given on 
U , and assume that there is an analytic involution system gi, . . . ,gs in the open 
subset U such that [/, gi] = for 1 < i < s . Then, given arbitrary analytic 
functions /ii, . . . on U , the vector field 




D'if{x)q{x) = 



is reducible by the common invariants of gi, . . . , g_ 



Proof. 



For any k one has 




due to [gk, f] = 0. The assertion follows from Theorem 12.71 
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Stating a converse to Proposition 12.111 is not a straightforward matter. 
While an involution system uniquely determines the analytic invariants, the in- 
variants do not determine a unique (finite) involution system. (Incidentally, this 
observation may be used to prove Frobenius' theorem; see for instance Hermann 
But the following statement holds. 

Proposition 2.12. Let gi,...,gs form an analytic involution system on an 
open set U such that the vector field f* is reducible by the common invariants of 
the Qi. Then locally, near any point of U where {gi, . . . ,gs) has maximal rank s, 
the module spanned by the gi has a basis cji = (^ijdj > with analytic functions 
Cij , 1 < ^, j < and there exist analytic functions /ij such that f '■= f* — Yl l^jdj 
satisfies [/, gi] = for all i . 

Proof. Use Lemma [23| with invariants Xi, . . . ,Xm and module basis e^+i, • • • , e. 
Then (with the same notation as in the proof of Lemma 12. 5p 



/ fi{xi, . . .,x„^) \ / \ 



fm (-^l ; • • • ; •^m) 





V 

provides the asserted decomposition. 







/ V*/ 



Corollary 2.13. In the special case s = 1, g = gi a vector field f* is reducible 
by the invariants of g if and only if there is a vector field f and a scalar function 
p such that [g, /] = and f* = f + pg- 

Proof. Using the notation of Proposition lTTTl we have g = erg and / = f*—fig; 
the assertion follows with p = pa. ■ 

Remark 2.14. (a) In the special case s = 1, g = gi, a direct proof of Propo- 
sition IKJB runs as follows: Assume [g, /*] = Pg for some function P, and make 
the ansatz f = f* — p g ■ Then [g, /] = if and only if Xg{p) + /3 = 0, and the 
latter has a solution near any non-stationary point of g (e.g. by the straightening 
theorem). Thus the function p can be determined explicitly whenever a straight- 
ening map for g is explicitly known. 

(b) As noted above, passing to a different module basis cji = Y'^ijdj "^"^^^ '^^^ 
change the reducibility conditions and properties but may affect other distinguished 
properties, like commutation of vector fields. For the case of one module generator 
g = gi this was discussed in Pucci and Saccomandi JT^ . o 



2.3. Orbital reducibility. For first-order ordinary differential equations one 
is not only interested in symmetries but more generally in orbital symmetries. 
Similarly, it is sensible to generalize from reduction to orbital reduction, as we will 
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do next. We first recall a characterization of orbital symmetry from [16], Lemma 



Proposition 2.15. Let the analytic differential equation x = f{x) be given on 
U . Assume that there is a rank m analytic involution system gi,...,gs in the 
open subset U, with n — m independent common invariants ipi, . . . ,ipn-m such 
that every common invariant can be expressed as an analytic function of the ipj . 
Assume that f is not an element of the module generated by the gj . Then 
is orbitally symmetric with respect to gi, . . . , gg if and only if there is an analytic 
function fi without zeros on U such that 



Correspondingly, we define: 

Definition 2.16. Equation ([1]) is orbitally reducible by the map if some 
equation x = ■ /(x) (/i analytic without zeros on an open subset U C t/) is 
reducible by ; in other words, instead of ([2]) the identity 



holds on U . 

Remark 2.17. (a) This is clearly a necessary and sufficient condition for solu- 
tion orbits of to be mapped by to solution orbits of x = h{x) . Equivalently, 
for every local first integral p of x = h[x) the pullback po^ is a first integral of 
X = f{x). 

(b) Via autonomization this definition extends to non- autonomous equations, and 
due to part (a) and Remark IKE this is the natural notion of a reducing map for 
non- autonomous equations; i.e., some map which sends solutions of a system to 
solutions of a system "in smaller dimension" . o 

Remark 2.18. While explicitly finding a reducing map to a one- dimensional 
equation is a hard problem, finding an orbital reducing map to any one- dimensional 
equation with nonzero right-hand side is locally trivial: For any h such that ho"^ 
is not identically zero, one may choose the factor p in a suitable way. But this 
observation is of little interest since it provides no information about solutions to 



Next we want to give a characterization of orbital reducibility in terms of 
Lie bracket properties. The critical argument in one direction of the proof is similar 
to [TB], Lemma 2.3. 

Theorem 2.19. Let the analytic differential equation x = f{x) be given on U , 
and let a rank m analytic involution system gi, . . . ,gs be given in the open subset 
U . Assume that there are n — m independent common invariants ipi, . . . ,ipn-m of 
the gi such that every common invariant can be expressed as an analytic function 



2.3. 



[f^f, 9j] = 0, 1 < j < s. 




(3) 
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of the ipj . Then ([T]) is orhitally reducible by the ipj if and only if there are analytic 
functions ai and Xij on U such that 

[9i, f] = tti/ + X] "^i^Sj, I <i< s (4) 

3 

Proof. In the setting of Corollary 12.91 we obtain 

[l^- f, 9i] = ^lij9j, l<i<s, 

3 

and therefore 

■ [/, 9i] = XgXfJ^) ■ f + ^ lij9j, l<i<s. 

3 

Division by fi yields (jl]). 

For the reverse direction assume that (j4]) holds on ?7. If Xf{(j)) = for 
every common invariant of the gj then / lies in the module generated by the 
Qj, due to the rank and independence conditions, and the bracket relation holds 
trivially. Otherwise, let ip be analytic such that 

X,,{^) = ... = XgM = 0, but Xj{^) ^ 0. 

Setting 

f* := — ^ — f 

one obtains 

[9uf*] = J2^*39^^ 1<^<^ 
j 

with analytic A*^ on U* := {x & U : Xf{ilj){x) ^ 0}. Indeed, for any i the 
commutation relation for gi and / implies 



and thus 
In view of 



X,^Xfi^)=a,Xfi^). 

1 



-f 



''XfW\ XfW 
the assertion follows with A*^- := Xij / Xf{ip) . u 

Remark 2.20. If the involution system consists just of g = gi then one obtains, 
as a particular case, the condition [g, f] = af+Xg. This leads to the X -symmetries 
of Muriel and Romero; see Q 0, 0, [7^ . Theorem 2.5 in JEj corresponds directly 
to the above Theorem in case s = 1; see also Section 2 of ITO^ . The theoretical 
framework was clarified by Morando JBj/ (Subsection 4-2 in particular) . These 
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observations presume the identification of higher-order ODEs with an equivalent 
first-order system, and restriction of g to infinitesimal point transformations. (See 
more on this in Section\^ below.) For this reason one could call [gi, . . . , gs) a 
system 0/ joint A -symmetries for the equation ([I]), /n |^ the name a- symmetries 
was chosen, to emphasize the focus on prolongations. o 

Remark 2.21. In the setting of Theorem \2.19\ the actual computation of an 
orbitally reduced system works as follows. There exists some fi such that the 
identities 

[f^- 9i] =^ltj9j, l<«<s, 
j 

hold, and by Corollary \2.9\ there exist analytic functions hi such that 

for 1 < i < n — m. While fi may not be explicitly known, one may turn to 
h* := hi/ hi for 1 < i < n ~ m , which can be expressed as functions of the ipi 
alone. Therefore the ipi define an orbit-preserving map from l^^j to y = h* (y) . o 

We finish this subsection with the counterparts to Proposition 12.111 ff . The 
proofs are immediate, in view of Theorem I2.19[ Proposition 12. 15[ Proposition |2]TT1 
Proposition 12.121 and Corollary 12.131 

Proposition 2.22. Let the analytic differential equation x = f{x) be given on 
U , and assume that there is an analytic involution system gi, . . . ,gs and analytic 
functions ai in the open subset U such that [/, gi] = Uif for 1 < i < s . Then, 
given arbitrary analytic functions fii, . . . , fig on U , the vector field 

f* ■■= f + '^f^t9i 
is orbitally reducible by the common invariants of gi, . . . , gs ■ 

Proposition 2.23. Let gi,...,gs form an analytic involution system on an 
open set U such that the vector field f* is orbitally reducible by the common 
invariants of the gi. Then locally near any point of U where {gi, . . . , gg) has 
maximal rank s, the module spanned by the gi has a basis gi = J2j^ij9j' '"'^^^ 
analytic functions aij , 1 < z < s , and there exist analytic functions fii such that 
f '■= f* — J2l^j9j satisfies [/, gi] = aif with suitable analytic ai, for 1 <i < s. 

Corollary 2.24. In the special case s = 1 , g = gi a vector field f* is orbitally 
reducible by the invariants of g if and only if there is a vector field f and scalar 
functions a, p such that [g, f] = af and f* = f + pg- 

2.4. Reduction by group invariants. In this subsection we consider some 
aspects of the "inverse problem" to determine all vector fields that are (orbitally) 
reducible by some given map. Following a guiding principle established in Olver 
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and Rosenau [13] (albeit in the context of partial differential equations), we con- 
sider reduction of ([T]) by invariants of some group, with the system itself not 
necessarily symmetric. Given a (local Lie) group of transformations, it may be 
difficult to determine all differential equations which are symmetric with respect 
to this group. By extension it it may be difficult to determine all differential equa- 
tions which are reducible by its invariants. But at least Propositions 12.11] and [2^22] 
provide a simple construction of reducible systems from symmetric ones, and for 
local one-parameter groups this construction yields all reducible systems. We ex- 
tend this result and show that for Lie algebras of compact and connected (linear) 
Lie groups, there is a method to construct all reducible vector fields on an open 
and dense subset. The underlying reason is the existence of a convenient repre- 
sentation (on an open-dense subset) for any vector field. This may be considered 
a consequence of the slice theorem (see e.g. Brocker and tom Dieck [1]), but we 
use a simple shortcut. The following results are an extension of [3], Example 2.5, 
where the Lemma and the first part of the Proposition were proven. The remaining 
statements are clear from the previous subsections. 

Lemma 2.25. Let G C GL{n, R) be a connected compact linear group, with 
invariant scalar product (-, ■), and denote by Q its Lie algebra. 
Denote by s the maximal orbit dimension of G and let Bi, . . . , Bg & Q be such 
that Biz, . . . ,BsZ are linearly independent in for some z, hence for all z in 
an open-dense subset. (In other words, s is the rank of the involution system 
generated by Q .) Then there exist algebraically independent polynomial invariants 
(Ti, . . . , (Jn~s; let their gradients qj be defined by 

D(Jj{x)y = {qj{x), y) . 

The Qj are G -symmetric (thus every transformation in G is a symmetry for 
x = qj{x) ), moreover 

d{x) := det {Bix, . . . , BgX, qi{x), qn-s{x)) 
is a nonzero polynomial, and every vector field f on U admits a representation 

fi^) = ^ai{x) BiX + ^l3j{x)qj{x) (5) 
which holds on U := {x G t/; 9{x) ^ 0}. 

Proposition 2.26. Let the hypotheses and notation be as in Lemma \2.25[ and 
let the vector field f be represented as in (|5]). 

(a) The vector field f is reducible by the invariants of G if and only if all (3j are 
G -invariant. This is equivalent to 

being G -symmetric. Thus every reducible system on a subset of U is obtained from 
a symmetric one via Proposition \2. 1 1[ 

(b) The vector field f is orbitally reducible by the invariants of G if and only if 
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there is an analytic u such that all j3j = v ■ j3j , with /3j G -invariant. This is 
equivalent to 

fix) ■= 

being orbitally Q -symmetric. Thus every orbitally reducible system on a subset of 
U is obtained from an orbitally symmetric one via Proposition \2. 2S[ 



Remark 2.27. (a) Note that Lemma \2. 25\ and Proposition \2. 26] provide a con- 
struction of reducible systems which does not require a priori knowledge of all 
symmetric systems. (Actually, finding all reducible systems here is less trouble- 
some than finding all symmetric systems.) 

(b ) The polynomial 6 is not necessarily G -invariant ( contrary to the statement in 
but generally 9 is the product of a G -invariant polynomial 6i and a polyno- 
mial $2 with the property that Biv, . . . , BgV are linearly dependent for all zeros v 
of 02 (m C";. o 

Examples, (a) Consider in the system 

X = f{x) = a{x) -xi + I3i{x) 





This representation corresponds to Proposition 12.26( a). with the 1 -dimensional 
subgroup G of 5*0(3,/?) generated hj g = {x2, —Xi,0) , on the set U defined by 
{xf + X2) ■ X3 7^ 0. The system is G-symmetric if, and only if, a and the Pi are 
functions of the generating invariants ai = xi x\ and (T2 = X3 alone. (This 
holds because the Lie algebra is abelian. Generally finding symmetric systems is 
a more involved matter.) The system is reducible by the invariants ai, o"2 of g if 
Pi and ^2 are functions of x\ + x\ and X3 alone , with a arbitrary. Assuming e.g. 
/32 7^ 0, the system is orbitally symmetric with respect to g if and only if a//32 
and Pi/ 132 are functions of ai and (T2 alone, and the system is orbitally reducible 
by the invariants of g if and only if P1/P2 is a function of ai and (T2 alone. 

(b) We consider G = 50(3, M). On the linear maps 

Bi{x) =1 xi I , B2{x) =1 I , B^ix) 

Xi 

span the Lie algebra Q , one has s = 2, and may choose Bi and B2 , since Bi z and 
B2Z are linearly independent for all 2 G f/ := {x; Xi 0} . The invariant algebra 
is generated by the polynomial a{x) := + + x|, and its gradient is equal to 
q{x) = 2x. One has 9{x) = 2xi(xf + Xg + x|) and on the set given by 6* 7^ any 
vector field can be written in the form 

/(x) = ai(x) ■ Bix + a2{x) ■ B2X + /3(x) ■ x. 

This vector field is reducible by a if and only if /3 is group-invariant, thus can be 
written as a function of a alone. There is no restriction on orbital reducibility by 
a, as was to be expected from Remark [2. 181 
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3. Higher order equations 

In this section we apply the results on first-order systems to ordinary differential 
equations of higher order, and thus gain a new perspective on the construction and 
reduction of equations admitting lambda symmetries. Recall the correspondence 
between higher-order equations and systems of first order: Given a (single) non- 
autonomous equation of order m + 1 > 1 , 

=p(t,a:,i;,...,x('"-^)), (6) 

its solutions correspond to solutions of the first-order system 

Xl = X2 

■_ (7) 

■^m— 1 •^m 

Therefore symmetries of the first-order system will send solutions of ([6]) to solutions 
of ([6]). In other words, orbital symmetries of the autonomous system 

1 

briefly x = P{x), (8) 

p{Xo,Xi,X2, . . .,Xm), 

will send solutions of IQ to solutions of (j6]) (up to familiar identifications). This 
point of view is proposed in the monograph by Stephani [13]. It may be worth 
noting (and has already been noted in special instances, e.g. by Nucci and Leach 
[TTj ) that conversely any first order system may locally be represented as a single 
higher order equation. 

Proposition 3.1. Let a nonautonomous first-order system 

z = q{t, z) on U CKx K'" 

be given, and let [t*, z*) such that q{t*, z*) ^ 0. Then there exist local coordinates 
t, Xl, . . . , Xm near (t*, z*) in which the system takes the form ([7]). 

Proof. With no loss of generality we have t* = and z* = 0. Consider the 
autonomized system 

We may assume that (5(0, 0) = (1, 1, 0, ... , O)*"". By a straightforward variant 
of the straightening theorem we may furthermore assume that locally Q{t, z) = 
(1, 1, 0, . . . , 0)*'' is constant. Now let 

m 

.7=2 ^ 



Xq 
Xl 
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and define 



a;,:=^=Xj-^(0), 1 < k < m. 
At z = the Jacobian matrix of (xi, . . . , XmY^ as a function of z equals 



/ 1 ■■■ \ 

1 ■■■ 

■■■ 1 

\ 1 ■■■ y 



and therefore we have a coordinate transformation (t, z) H- {t, x) . By design one 
has Xj = Xj+i for 1 < j < m — 1 , and this proves the assertion. H 

Remark 3.2. (a) The proof is constructive, to some point, even if Q is given 
in general form (only assuming Q{0, 0) 7^ Take any function cf) such that 
(J), Xq{(I)), . . . , Xq~^ {(p) are functionally independent, and choose the new vari- 
ables accordingly. The above proof amounted to verifying the existence of such a 
function, and also makes clear that "almost every" function will satisfy this prop- 
erty. More precisely, in the local ring of analytic functions those which do not 
satisfy the independence property form a subset of positive (finite) codimension. 

(b) The exceptional case of functionally dependent (p, Xq^cj)), . . . , Xq~^ ((p) pro- 
vides reduction in a direct manner. Let i be maximal such that Xi := 0, X2 : = 
Xq(0), ...,X£ := are independent. Then Xq^cf)) may be expressed as a 
function of Xi, . . . ,xe, whence one has reduction to an equation of order i+1 (or 
the associated system). 

(c) This Proposition opens, in principle, a possibility to determine symmetries of 
first-order systems: Rewrite the system as a higher-order equation, and determine 
the point symmetries of the latter, which amounts to a fully algorithmic procedure. 
(Nucci and Leach [77] / noted and used a variant of such an approach.) The draw- 
back is, of course, that one will only find trivial symmetries in general. It could be 
interesting to explore the possibility of a systematic approach, which would include 
the question how to choose a suitable function cf). o 



3.1. Review: Prolongation in simple cases. 

In most monographs on symmetries, such as Olver [12], or Krasil'shchik 
and Vinogradov [S], the Lie point symmetries of higher order equations (E]) are 
determined via the general procedure for prolongations of vector fields to jet spaces. 
The following shortcut works and is quite useful in our scenario; it is essentially 
taken from Stephani [15], Ch. I, Section 3.5. 

We are interested in vector fields g on U such that [g,q\ = IJ' ■ q for 
some n . By geometric motivation one frequently considers only infinitesimal point 
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symmetries of the higher order equation, which imphes the projectabihty property 



g2{Xo,Xi,X2) 

Qm—l iS^Qi ■^li • • • ; -^m— l) 
y gm{xQ,Xi, . . .,Xm) J 



(9) 



for the infinitesimal symmetry of the associated system. Step-by-step evaluation 
of the condition [g^q] = fJ' ■ q now yields 



- ^ofi-o - X2 ■ digo 
g2 - dogi - X2 ■ digi 
gs - dog2 - X2 ■ dig2 - X3d2g3 



/iX2 
fJ,X3 



gr, 



dog 



m— 1 



X2 ■ dig, 



m — 1 



Xmdmgm—1 



llXr, 



Thus one may successively compute yU, g2,...,gm from go and gi and their 
derivatives. Evaluating the last entry of the Lie bracket, which has not been 
written down here, provides an overdetermined, algorithmically accessible, system 
of partial differential equations for g^ and gi, and thus in effect the symmetry 
conditions. (This approach is of course equivalent to the usual prolongation 
procedure for point symmetries of (|6]).) 

Generalizing to orbital reducibility, but keeping the geometric restriction 
one deals with the lambda symmetries first considered by Muriel and Romero 
[7]. The condition is 

[g,q] = fi - q + X- g. 

Initially Muriel and Romero [7] (for geometric reasons) require A to be a function 
(possibly a priori unknown) of xq, Xi, X2 only. One obtains the " A -prolongation 
formulas" : 



- dogo - X2 ■ digo 
g2 - dogi - X2 ■ digi 
gs - dog2 - X2 ■ dig2 



X3d2g3 



fi + Xgo 
11x2 + Xgi 
fixs + \g2 



gr, 



dog 



m— 1 



X2 ■ dig 



m— 1 



• • • XjYiOjyigm — 1 



fiXm. + \g 



m— 1 



In this setting. A, 5^0 and gi successively determine /i, g2, . ■ ■ ,gm, and again the 
last entry of the Lie bracket identity will provide compatibility conditions, as was 
noted by Muriel and Romero [7]. The determination of A -symmetries is not a 
completely algorithmic procedure, and the artwork in this approach is to suitably 
determine (e.g. by educated guesses) A such that prolongation and evaluation 
yields nontrivial results. Muriel and Romero's restriction imposed on A is of 
importance for the construction of higher order differential invariants further on 
in [7]. One should also note that these authors relaxed the geometric restrictions 
on A in a subsequent paper [5], thus making another step towards general orbital 
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reduction with respect to a single vector field, albeit in the setting of a higher 
order equation. (See also Remark [2.201 ) 

In the following discussion and construction of reducible higher-order dif- 
ferential equations we emphasize the correspondence to first order systems. 

3.2. Constructing reducible higher-order equations. We propose a con- 
struction of reducible higher-order equations, thus extending the work by Muriel 
and Romero [3 El El [ID] on lambda symmetries. We use Theorem 12 . 1 9 1 and Proposi- 
tion |222l based on the correspondence to first-order systems established in Propo- 
sition [3]T1 starting from a system with known orbital symmetries. Thus, on the one 
hand, as in Muriel and Romero [9] the "lambdas" may depend on all variables, and 
moreover we do not necessarily restrict attention to point symmetries of higher or- 
der equations. On the other hand, we extend the framework of Muriel and Romero 
by considering more than one infinitesimal orbital symmetry. As noted earlier, we 
therefore work in the most general setting for (orbital) reducibility of the asso- 
ciated first order system. In contrast to Muriel and Romero [71 El [10] our focus 
is on constructing reducible equations, rather than detecting reducibility in given 
equations. We emphasize that the following should be seen only as a first step 
towards a systematic construction of reducible equations with prescribed (joint) 
A -symmetries. 

Consider a single equation (|6|) of order m + 1 and rewrite it as the autono- 
mized system ([H]): 

±0 =1 

Xl = X2 

Xm P('^0; "^l; • • • ; •^m} 

Assume that for k = 1, . . . ,r the system admits infinitesimal symmetries 



9 



/^f(xo,x„.., 



\ gi!;\xo,Xi, 



which form an involution system. Then for any system of scalar functions u^'^^ the 
equation 

f i + E.^^^)^?f(a:) \ 

X2 + Ek''^''^9['\x) 



X = H(x) 



is orbitally reducible by the common invariants of the g^^^ . Now consider the 
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orbit ally equivalent system 



X = H(x) 



1 + E^(^)^7f I 



-H{x) 



X) 



and introduce new coordinates 



t = xq, yi = xi, 1/2 = hi{x), 



( 1 \ 

hx{x) 



which is always possible by Proposition 13.11 In general one will thus obtain an 
equation of order m + 1 ioT y = yi. (The exceptional case when the functions 
Xi, Xh{xi), . . . , X^~^(xi) are not independent implies reducibility via Remark 
13.21 ) By construction, the equation x = H{x) is orbitally reducible by the 
common invariants of the g^''^ , and in the generic case this holds true for the 
system in new coordinates t,yi, . . . ,ym (with the invariants also written in new 
coordinates), and for the corresponding equation of order m + 1. Thus we have 
constructed reducible higher order equations from symmetric ones. Note that if 
the coordinate transformation is given by a; = "^{y) then the "joint- A symmetries" 
in new coordinates are given by 



-'g^'\^{y)). 



Remark 3.3. The special case when r = 1 and the geometric restrictions 
on g = gi ( and A ) hold is, naturally, of particular interest. It is possible to 
explicitly construct all X -symmetric higher-order equations which are reducible 
by the invariants of g, assuming the latter are known. First, one knows that 
all symmetric higher- order equations are given by functions of the differential 
invariants of g (with appropriate identifications; see e.g. Olver fW^ . Ch. 2), and 
second, by Corollary \2.24\ one can construct all reducible equations via the procedure 
outlined above. Here one should assume v to be a function of Xq, Xi and X2 only, 
so that the same holds for A . o 



Remark 3.4. In the special case r = 1 it is also of interest to identify the 
"X " emerging from this procedure. Thus start with vector fields f and g such that 
[Q-i /] = ct/ with some scalar function a. Given a scalar function v , form 



H{x) = f{x) + ug{x); H{x) 



Straightforward computations show that 



-H(x). 



aH + {Xg{u) - au) g; [g,H] 



i^go[x) 



:...)H 



and the coefficient of g in the second identity (rewritten in new coordinates via 
X = '^{y)) is the scalar function X as introduced in Q/. o 
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3.3. Examples of order two. Since orbital reduction to dimension one is of 
little interest (recall Remark I2.18p . "joint A" is of little interest here. Therefore 
we remain mostly within the framework of Muriel and Romero [3 [9], considering 
the inverse problem of finding differential equations with prescribed reduction. 
Rewrite a single second order equation as an autonomous system: 

±0 = 1 

Xi = X2 

±2 = p(Xo,Xi,X2) 

Assume that 

go{xo,xi,X2) 
gi{xo,xi,X2) 
5(2(3:0,2:1, xa) 

is an infinitesimal orbital symmetry for this equation. Then for any scalar function 
u the equation 

/ 1 + ugo{x) 
X = H{x) := X2 + ygiix) 
\ p{x) + ug2{x) 

is orbitally reducible by the invariants of g, and every orbitally reducible system 
is obtained from an orbitally symmetric one in this way, due to Proposition 12.221 
and Corollary 12.241 The orbitally equivalent system 

X = Hix) := i-^Hix) 

l + ugo{x) 

remains orbitally reducible by the invariants of g. Now, unless hi depends on xq 
and Xi alone, there is a local coordinate change 




2/0 = Xq 

yi = xi 

2/2 = Xnivi) = hi{x) 
and the system in new coordinates 



y = H*{y) = 




is orbitally reducible by the invariants of g, expressed in new coordinates. Thus 
we have obtained a reducible second-order equation from a symmetric one, and 
every equation which is orbitally reducible by the invariants of g is obtained in 
this way. Generally, the method is not completely constructive, since an explicit 
computation of p* requires an explicit inverse to the coordinate transformation. 
Therefore we (have to) make special choices of functions in the concrete examples 
below. 
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Example 1. Consider the system 

Xq = 1 

Xi = X2 (10) 

Here 7 is an arbitrary analytic function of one variable. This system admits the 
infinitesimal orbital symmetry 




(which corresponds to a Lie point symmetry of the associated second-order equa- 
tion). An independent set of invariants of g is given by 

CTi = Xi/Xo (72 = X2. 

Therefore, given any scalar function u the equation 

(1 + l^Xo 
X2 + lyxi 
1{X2)/X1 

is orbitally reducible by these invariants, and indeed one has: 

Xo&i := Xo ■ X^((Ti) = (cTa - c^i) Xoa2 := Xq ■ ^^^(aa) = 7(o-2)M 
i.e. Definition 12. 161 applies with fi = xq. The orbitally equivalent system 



X = Hix) = ^H(x) 




1 

X2 + UXI 
l + UXQ 



l + iygo{x) \ j^j^^ ] \ 7(^2) 

remains orbitally reducible by the invariants of g . The coordinate change is now 



Vo = Xo 
yi = xi 

nr^i — 



y2 = hix)-^^^ 



(i) In the particular case of constant u, the system in new coordinates is 

y = H*{y) =\ y, \:=\ y2 | (11) 




where 

X2 = ^2(1 + i^yo) - i^yi- 
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This system is orbitally reducible by the invariants of g , expressed in new coordi- 
nates, and we note that A = u/i^l + u) according to Remark 13.41 The invariants 
are 

5'i = 2/1/2/0 5^2 = 2/2(1 + i^2/o) - ^^2/1 
and the orbital reduction is obtained with Jl = yo{l + uyo) . In detail 

Jl-Xn'-ioi) = (5'2^-^i) 
Jl-XH*{a2) = 7(5'2)/5i 

Choose now for instance 7 = Then Equation ffTOj) gives the ODE xx = x^ 
which is easily solvable, whereas the ODE resulting from f lTTj) (with xq = yo = t 
and yi = y) is 

yy{l + vtf = {y + vtij - vyf (13) 

which seems to be not solvable by standard methods. However, from ( 1T2|) with 
7 = ^2 deduces 

5^2 ~ 

— log (J2 = const 

0"! 

which expresses a first integral for system flTT]) . On the other hand, a first integral 
for this system corresponds to a first integral for the resulting second order ODE 
(|T3|) . Indeed, it can be checked that 

— (y + lyty — uy) — \og(y + uty — vy) = const 
2/ 

is satisfied. 

(ii) In the particular case that u = X2 one has 

y2 = X2[ — , X2 = — =: 2/ , A = — ■ -. 

VI + 0:2X0/ 1 + 2/1-2/02/2 1 + 2/1 + 2/2(1-2/0) 

The invariants, expressed in new coordinates, are now yi/yo and 0(2/). From the 
third entry of H* one sees that 

7(0) ■ 2/2 , , wi ■ + 7(0) ■ yo 

= T~ + ■ 2/2 - ■ , , — ^ 

2/1-0 (1 + ■ yoV 

As a particular example with 7 = we obtain the second-order equation 



y ^l-y-ty (1 + 2/)^'^' 



• 2 



which is reducible to a first order equation via the invariants of g . 

Example 2. The previous example still remains in the classical setting of lambda 
symmetries, since we started with a point symmetry. For an example in a more 
general setting, start with the simple system 

Xq = 1 

±1 = X2 (14) 
X2 = 
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which corresponds to i; = and admits the infinitesimal symmetry 

9 

Given any scalar function u the equation 

X = H(x) 




1 + 1^X2 




is orbitally reducible by the invariants of g . An independent set of invariants of g 
is given by 

O'l = Xq- X2 02 = X2/X1. 

Let us consider the special case with z/ = xi . The coordinate transformation then 
is given by 

xj + X2 y2- yl 

yo = xo, yi = xi, y2 = — , ^2 = - 

1 + X1X2 1 - yiy2 



and we have 

A 

In the X -coordinates we find 



1/1(1-^11/2) 

i-yf ' 



Xh{(Ti) = 1/(1 + 0:1X2) 
Xh{(^2) = -ai/(l +X1X2) 

and thus we have orbital reduction of x = H{x) to the system 

ii = 1 

Z2 = —Z2- 



(15) 



We have a reduced (autonomous, in this particular case) one-dimensional equation 
dz2/dzi = —z^. 

Going to new coordinates the system x = H{x) is equivalent (by straight- 
forward computation) to the second-order equation 



{yy -y^) n , .^2 , y-y^ + '^yy -f- y^y^ 



yy) + 



which therefore is reducible to Equation f[T5|) by the invariants ai and 5^2 (ex- 
pressed in new coordinates). Note that x = H{x) is autonomous and therefore 
admits time translation as a symmetry, which is refiected in the one-dimensional 
orbitally reduced equation also being autonomous. (We chose this system for the 
sake of brevity, but note that the reduction to (IT^ is not a symmetry reduction.) 

3.4. An example of order three. 
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The essential point of this example is to illustrate nontrivial "joint lambda" 
symmetries. We start with the simple system 



Xo 

is 



1 




corresponding to the third-order equation x^"^^ 
infinitesimal orbital symmetries 



(16) 



0. This system admits the 



/ ^0 \ 

Xi 



V j 



.(2) 



/ \ 

Xx 
X2 
\X3 / 



[9^'\ 9^'^ = 0, [g^'K f] = /, [9^'\ f] = 0. 



Note that g^^^ and g^'^'^ both are point symmetries. A set of independent common 
invariants of g^^^ and g^'^^ is given by 

■01 = X0X2/X1, ^2 = xlx-i/xi. 

We construct an equation that is reducible by the common invariants of g^^^ and 
g^'^K Let i/^^\x) = Xi and z/(^^(x) = 1/xi, thus 



H(x) 



( 1 + XqXi \ 
X2 + + 1 



and 



E{x\ 



1 

3:2+3:1+1 

X1X3+X2 



\ 



2:1(1+2:0x1) 

. (l-2:f)a:3 . 
\ a;i(l+xoxi) / 



Let us look at the reduction first. One finds 



^5(^1) 



and therefore ^^l and '?/^2 provide an orbital reduction of H (as well as of H) to 
the autonomous two-dimensional system 



Z2 



Zl - Z( + Z2 
2Z2 — Z1Z2, 



(17) 



which may be rewritten as a non-autonomous first order equation: 

2^2 ~ Z1Z2 



dZ2 

dzi 



Zl — zf — 2Z2 
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The example was primarily chosen to obtain an explicitly invertible coordinate 
transformation towards the third-order equation. Since some expressions are some- 
what unwieldy, we will write them down only in an abbreviated version. Passing 
to new coordinates, we set 

Vo = xo, yi = xi, y2 = XHiVi) = — h 



1 + yovi 1 + yoyi 
which yields 

X2 = 1/2(1 + 2/02/1) - (1 + 2/i)- 
Using the "hybrid" expressions 

1 + 2/02/1 2/1(1 + 2/02/1) 2/1(1 + 2/02/1) 
routine calculations provide 

-yl - yoym + 2?/^?/2 -yi + yoyfy2 + (1 - - ymy2)x2 + yixi. 



2/3 



2/1(1 + 2/02/1) 



which (being linear in 0:3) can easily be solved for as a function of yo, ■ ■ ■ ,y3- 
Taking the Lie derivative Xniys) and making the usual identifications, one obtains 
a (lenghty) third order equation for y = yi, which can be reduced to (1171) by the 
invariants tpi and ip2 (expressed in the y-coordinates). 

The example shows that the construction of nontrivial reducible higher 
order equations via prescribed "joint lambda" symmetries is feasible. But it also 
illustrates that work remains to be done towards a systematic approach. 
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